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We propose using matrix-isolated paramagnetic diatomic molecules to search for the electric 
dipole moment of electron (eEDM). As was suggested by Shapiro in 1968, the eEDM leads to a 
magnetization of a sample in the external electric field. In a typical condensed matter experiment, 
the effective field on the unpaired electron is of the same order of magnitude as the laboratory field, 
typically about lO^V/cm. We exploit the fact that the efi^ective electric field inside heavy polar 
molecules is in the order of 10^'^V/cm. This leads to a huge enhancement of the Shapiro effect. 
Statistical sensitivity of the proposed experiment may allow one to improve the current limit on 
eEDM by three orders of magnitude in few hours accumulation time. 

PACS numbers: 11.30.Er, 32.80.Ys 



The searches for the elusive electric dipole moment of 
electron (eEDM) are motivated by the fact that the ex- 
istence of a permanent EDM of a particle violates both 
parity (P) and time-reversal (T) symmetries. Due to 
the compelling arguments of the CPT theorem, the T- 
violation implies CP-violation, a subject of great interest 
in the physics of fundamental interactions Current 
experimental limit on eEDM is close to the predictions 
of many extensions to the Standard Model of elementary 
particles, such as "naive" supersymmetry (SUSY) pUj]. 
Other SUSY extensions yield eEDM a few orders of mag- 
nitude below the present limit. Here we propose an 
eEDM search that may constrain eEDM at that impor- 
tant level. Our proposal relies on the fact that the ther- 
modynamically averaged eEDM (and thus the electron's 
magnetic moment aligned with eEDM) is oriented along 
the electric field. We propose to employ polarized molec- 
ular radicals frozen in a rare-gas matrix and measure the 
eEDM-induced magnetic field generated by the sample. 
Conservative estimates project that the present limit on 
eEDM can be improved by several orders of magnitude. 

The present limit on eEDM, 

141 < 1-6 X 10"2^e-cm, (1) 

is derived from a high-precision measurement |3| with a 
beam of Tl atoms. In such experiments one spectroscop- 
ically searches for a tiny eEDM-induced splitting of the 
magnetic sublevels of an atom in an externally applied 
electric field. 

New atomic eEDM experiments plan to use optical 
trapping There are two other major trends aimed 
at improving the experimental sensitivity to eEDM: (i) 
employing molecules instead of atoms in spectroscopic 
experiments 0, 0^_@L^'^'i (ii) non-spectroscopic sohd 
state experiments |2lUIL Hj- Here we propose to merge 
these two trends by searching for eEDM with molecules 
trapped in a cold matrix of rare-gas atoms (see Fig.^l. 



We argue that this scheme combines advantages of both 
techniques. Indeed, the eEDM effects in molecules are 
markedly amplified because of the strong internal molec- 
ular electric field ^3|> much larger than attainable lab- 
oratory fields. In the present solid-state schemes the 
atomic enhancement of the external electric field for ions 
of a solid is of the order of unity ^3 ■ By using matrix- 
isolated diatomic radicals, one can gain up to six orders 
of magnitude in the effective electric field. At the same 
time one retains a great statistical sensitivity of the solid- 
state searches. We show that this particular combina- 
tion seems to drastically improve sensitivity of the eEDM 
search. 

Let us review important aspects of the non- 
spectroscopic solid-state search for the eEDM. Intro- 
duced by Shapiro this scheme exploits the link be- 
tween EDM of the electron and it's spin, d = deCr, 
and therefore it's magnetic moment, fie ~ —fJ'BC — 
—fiBd/de- In an external E-field, because of the cou- 
pling of the eEDM to the E-field, thermal populations of 
the spin-up and spin-down states slightly differ, leading 
to the magnetization of the sample. By measuring the 
generated magnetic field one derives constraints on the 
eEDM. A proof-of-concept experiment has been carried 
out in 1978 by Vasiliev and Kolycheva 0|. At that time 
the solid-state experiment appeared to be less sensitive to 
the eEDM than the spectroscopic beam experiments. It 
is only very recently that the advances in magnetometery 
(see [l5| and references therein) have revived an interest 
to the solid-state eEDM searches H ES Elli3- Alter- 
natively, one can look for a voltage induced in a sample 
in external magnetic field 

We focus on molecular radicals (i.e., molecules with un- 
paired spin) in the ground state. Consider a sample 
of radicals in thermodynamic equilibrium at temperature 
T. Because of the eEDM coupling to internal molec- 
ular E-field, spin substates in a molecule have slightly 
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FIG. 1: Scheme of searching for EDM of electron with di- 
atomic radicals embedded in a matrix of rare-gas atoms. A 
polarizing electric field E is applied to the matrix. As a re- 
sult, molecular CP-violating magnetic moments fjP^ become 
oriented and generate ultraweak magnetic field B'^^ . By mea- 
suring B'^^ one places constraints on eEDM. 



different energies. This mechanism leads to a thermo- 
dynamically averaged CP-violating (P,T-odd) magnetic 
moment per molecule (/i*"^^) ~ jiB deEcS / {ksT), where 
EeS is the large molecular effective electric field act- 
ing on the EDM of the unpaired electron. EeS grows 
(X with the nuclear charge Z of the heavier molec- 
ular constituent Il9| and one would choose to work 
with heavy radicals. Such molecules as BaF, YbF, HgF, 
and HgH belong to this broad category. We found that 
mercury hydride (HgH) has parameters most suitable for 
the proposed search, see Table For the HgH molecule 
£;off ~ 8 X 10i"V/cm 113 and its ESR spectrum in Ar 
matrix has been studied in 1211. 



TABLE I: Parameters of several heavy molecules with the 
ground state '^Si /o. Molecular dipole moments D were mea- 
sured in Refs. |2^. Polarization (riz), the maximal 
number density rimax, and the accumulation time tacc required 
to reach S/N=l for the current limit on eEDM Q. These pa- 
rameters are calculated with the help of Eqs. O, ©, and 11211 
for E = IQ kV/cm, T = IK, and sample volume 0.1 cm^. 



Molecule 


IT' a 
-C'eff 


D 




^max 


^acc 




(10''^) 


(D) 




(10^°^^ 


(ms) 


BaF 


8 


3.17 


0.13 


0.03 


300 


YbF 


26 


3.91 


0.16 


0.02 


30 


HgH 


79 


0.47 


0.02 


1.5 
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"The effective electric field for BaF and YbF was calculated in 
HMUBl' For HgH we rescale results from |2(3| . 

For diatomics, the moment {iJp^) is directed along 
the molecular axis. For a randomly oriented sample, 
however, the net magnetization would vanish. When 
an external E-field is applied, it couples to the tradi- 
tional molecular electric-dipole moment D and orients 
the molecules. Taking into account molecular polariza- 



tion, the CP-moment can be expressed as 



^J.B■ 



(2) 



where (n^) is the average projection of the molecu- 
lar axis onto the E-field (the field is directed along z- 
axis). Now the sample acquires a macroscopic magneti- 
zation. This magnetization generates an ultraweak mag- 
netic field proportional to eEDM 



(3) 



where n is the molecular number density and 7 is a 
geometry-dependent factor. For example, near the cen- 
ter of a disk-shaped sample of radius R and thickness L, 
7 = L/2R. 

Orientation of B-field © is linked to that of the ap- 
plied E-field through (riz). Such a Hnk is forbidden in the 
traditional electrodynamics. Its very presence is a man- 
ifestation of the parity and time-reversal violation. By 
measuring B'~^^ one constrains eEDM via Eqs. (jSJ and 
0. It is apparent that maximizing n is beneficial. How- 
ever, bringing radicals together is problematic — they 
react chemically. Here is where the matrix isolation tech- 
nique jl^l becomes key. In this well-established method, 
the molecules are co-deposited with rare-gas atoms or 
other species onto a cold (T ~ IK) substrate and be- 
come trapped in the matrix. Small trapped molecules 
exhibit properties similar to those for free molecules and 
a variety of studies, including determination of hyperfine- 
structure constants has been carried out. 

There is an upper limit on the density of trapped 
molecules; to avoid spin alignment in the subsystem 
of guest molecules one requires that thermal agitations 
are stronger than dipole-dipole interactions between the 
molecules. We can estimate the maximum density as: 



3 ksT 



An 



(4) 



A particular advantage of HgH is that its dipole moment 
is relatively small, D = 0.47Debye [lil and at T = 1 K, 
the density rimax w 1.5 x 10^°cm~^. 

Estimate (0J agrees with experimental observations 
that 1:100 guest to host ratio is possible. According to 
[2^ the realistic matrix thickness and area are L = 1 mm 
and S = 1 cm^. That corresponds to 7 « 0.1 in Q. Re- 
cently developed low density plasma beam source |2^ 
produces permanent beam of heavy radicals with inten- 
sity ^ 10^^ mol/sterad/s. Placing 1 cm^ target at 20 cm 
from the source, one can accumulate necessary number 
of radicals, i.e. 10^^, in 1 hour. 

How are the relevant molecular properties modified by 
the matrix environment? A free non-rotating molecule 
may be described by the electronic wave function |0), 
with O = ±1/2 characterizing projection of spin onto 
molecular axis. The time-reversal operation T converts 
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ri-states into each other: jfJ) ^ | — fi). In the matrix, a 
molecule can be considered as an individual entity per- 
turbed by the host atoms. The local symmetry of the 
perturbing fields depends on the position of the molecule 
in the matrix. Independent of the spatial symmetry 
the time-reversal symmetry remains. According to the 
Kramers' theorem, in the absence of magnetic fields, all 
levels of diatomics with half-integer spin remain two-fold 
degenerate for any possible electric field. 

EDM interaction operates at short distances near the 
heavier nucleus. Expanding the electronic wavefunction 
in partial waves we notice that contribution to the eEDM 
signal of total angular momenta beyond Si/2- and Pi/2~ 
waves are strongly suppressed because of the growing 
centrifugal barrier and properties of the eEDM • The 
truncated wave function has the Coc,v symmetry and f2 
still remains a good quantum number for the degener- 
ate states of matrix-isolated radicals. Within this ap- 
proximation, the effective molecular Hamiltonian in the 
external field E reads 



D E* 



(5) 



where E* is microscopic E-field; for small fields E* = 
E/e. We used Hcs to arrive at Eq. 

Using the estimate (0) with the present limit on 
eEDM we obtain for the thermally- induced CP-odd 
magnetic moment of HgH molecule trapped at T = 1 K 



(/zCP(HgH))<1.4xlO-i2(n,)MB. 



(6) 



It is instructive to compare this value to the permanent 
molecular CP-violating magnetic moment introduced by 
us in Ref. [T^ . This moment arises due to a magneti- 
zation of the molecule by its own electric field (irrespec- 
tive of the temperature). The largest /i^^ for diamag- 
netic molecules was found for BiF for which /x'-'^ < 3 x 
10~^'^ {nz)fiB, much smaller than the thermally-induced 
CP-odd moment ©. Therefore here we may neglect the 
permanent . 

An important parameter entering (/in^oi) is the degree 
of molecular polarization (n^) in the external E-field. 
Free diatomic molecules can be easily polarized by the 
laboratory fields ~ lO^V/cm, but there is a paucity of 
data on polarizing matrix-isolated molecules [2jJ. Cer- 
tainly, the rotational dynamics of the guest molecule is 
strongly affected by the matrix cage. The molecular axis 
evolves in a complex multi-valley potential, subject to 
the symmetry imposed on the molecules by the matrix 
cage. Depending on the barrier height between different 
spatially oriented valleys, the guest molecule may either 
execute hindered rotation or librations about the valley 
minima. Ref. |32l| reports evidence for hindered rota- 
tion of HXeBr and Ref. |3^ suggests that other hydrides 
can rotate. Note also, that for Ar matrix the cell size is 
4.5A, while internuclear distance for HgH is only 1.7A. 
That gives us a confidence that the HgH radical can be 
polarized by the external electric field. 



We will distinguish between two limiting cases of 
molecular polarization: strong and weak fields. In the 
former limit {nS) ~ 1, and in the latter, 



, , 1 fDE*n, 
{riz) = 2^n^exp 



DE* 



(7) 



For isotropic orientational distribution, characteristic for 



the polycrystalline matrixes. 



1/3, and we get 



CP \ 
mol/ 



1 DE* E,sde 
3 fcsT ksT 



(8) 



The dielectric constant of the rare-gas matrix is close to 
unity, but addition of polar molecules results in 



1 -|- 47rna = 1 + 47m- 



1 + — n , (9 

3 keT' ^ ' 



where a is molecular polarizability. For maximum den- 
sity Q, e « 2 and E* « E/2. 

The parameter differentiating the weak- and the high- 
field regimes is the ratio DE* /ksT. For HgH trapped at 
T = 1 K, the transition occurs at E* « 100 kV/cm. The 
breakdown fields for the rare-gas matrices are unknown, 
we only notice that for liquid Xe it is 400 kV/cm so 
that both weak- and high-field regimes may be possibly 
realized. The moderate E = 10 kV/cm field corresponds 
to (n^) w 0.02. 

Finally, we proceed to evaluating the sensitivity of the 
proposed eEDM search. There are two crucial criteria 
to consider: weakest measurable B-field and signal-to- 
noise ratio. Presently, the most sensitive measurement 
of magnetic fields has been carried out by the Princeton 
group (see |3| a-nd references therein). This group has 
reached the sensitivity level of 5.4 x 10~^^G/VHz. A 
projected experimental sensitivity of 3 x 10^^^G/\/Hz 
is published in 0. We find that for (n^) ~ 1 and for 
7 = 0.1 the present eEDM limit may be recovered within 
integration time of i = 5 s for the demonstrated sensitiv- 
ity and within 10^^ s for the projected sensitivity. Alter- 
natively, during a week- long measurement, the present 
eEDM limit may be improved by 3 x 10^ for the demon- 
strated and by 6 X 10^ for the projected B-field sensitivity. 
These values are reduced by a factor of 50 for a moderate 
10 kV/cm polarizing field. 

In addition to limitations imposed by the weakest mea- 
surable B-field one must also consider signal-to-noise ra- 
tio 34]. As we pointed out above, the thermally-induced 
(Mmoi) radicals is many orders larger than perm anent 
A*moi of diamagnetic molecules discussed in jl3|- The 
magnetic noise from paramagnetic radicals is also much 
higher as they have traditional magnetic moments asso- 
ciated with unpaired electron spin. 



(10) 



These moments lead to random magnetization of the 
sample and generate a fluctuating B-field. Unlike _B*-^^, 
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this field is not correlated with the direction of the ex- 
ternal E-field and it is the main source of the noise. In 
our case, the signal-to-noise ratio is 

S/N ^'3^^^^\lAf^ , (11) 
Mb V T 

where Af is the number of molecules, t is the observa- 
tion time, and r is the correlation time for the random 
thermal magnetization. Factor 3 at the right hand side 
appears because of the averaging of the magnetic moment 
(ll()|l over orientations of the molecular axis n. 

For a strong spin-rotation coupling, as in the case 
of HgH, r is determined by interaction of molecular 
axis with environment. One of such mechanisms is 
the dipolar interaction between guest radicals, so that 
T ^ h/{D'^n) = 47r?i/(3fcsr) for the optimal density (gj. 
For the weak-field limit JHJ we get the final expression 
for S/N: 

s/N = ^^^VvTJn, (12) 

on kb^ 

where V is the sample volume. This equation is used in 
Table U to estimate accumulation time needed to repro- 
duce the current limit ^ . For HgH molecule we find that 
for a volume of 0.1 cm^ and strong polarizing field, the 
present eEDM limit may be recovered within t = 10~^s 
(3 ms for the field 10 kV/cm). By integrating the signal 
for one week, the present eEDM limit may be improved 
by a factor of 2 x 10^. Note that these estimates are close 
to the estimates based on the projected sensitivity to the 
weak magnetic fields 0. 

To summarize, our proposed eEDM search combines 
advantages of the strong intermolecular field with a high 
attainable number density of molecules embedded in a 
matrix of rare-gas atoms. We argue that our proposal 
has a potential of improving the present eEDM limit by 
several orders of magnitude. That will allow constraining 
the "new physics" beyond the Standard Model at an im- 
portant new level and, in particular, testing predictions 
of competing SUSY models. 
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